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RESULTS  OF  DETAILED  FLOW  FIELD 
AND  RATE  CHEMISTRY  CALCULATIONS 


ON  AN  AERO BALLIST IC  PELLET 


By  Arthur  Fields 


I.  INTRODUCTION 


L 


The  method  of  characteristics  as  applied  to  a  real  gas  (assuming 


frozen  composition  along  streamlines)  was  used  to  compute  the  detailed 
supersonic  flow  field  about  the  pellet  geometry  shown  in  Figure  1.  The 
free  stream  velocity  was  22,  100  feet  per  second  while  the  free  stream 
pressure  and  temperature  were  1  mm  Hg  and  300°K  respectively. 

The  calculation  of  the  flow  field  properties  in  the  subsonic  region 

/■ 

was  performed  by  applying  the  technique  given  in  Reference  1  to  a  real 
gas.  The  air  was  assumed  to  be  in  thermal  and  chemical  equilibrium. 
The  results  of  this  calculation  served  as  input  to  the  frozen  character¬ 
istics  program  which  handles  the  computation  of  the  entire  supersonic 
flow  field  region  up  to  the  near  wake  shock.  Details  of  this  frozen 
characteristics  program  may  be  found  in  Reference  2. 

While  both  the  subsonic  and  supersonic  analyses  have  been  pro¬ 
grammed  for  the  IBM  704  computer  and  are  completely  automatic,  the 
calculation  of  the  wake  shock  was  carried  out  by  use  of  a  Bendix  G-15 
computer.  This  calculation  requires  monitoring  because  of  the  limited 
storage  capacity  of  the  computer.  The  flow  field  downstream  of  the 
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wake  shock  is  calculated  by  the  method  of  characteristics  utilizing  the  IBM 
704. 

Thus,  having  completed  the  calculation  of  the  detailed  flow  field 
about  the  pellet,  we  then  proceeded  to  study  the  effects  of  rate  chomietry. 
Because  of  the  complexity  of  solving  for  the  general  flow  field  properties 
of  a  gas  in  non-equilibrium  molecular  dissociation,  a  one  dimensional 
approach  to  the  solution  of  this  problem  was  used,  ( Roferonce  4. ) 

The  solution  of  the  detailed  frozen  flow  field  enables  us  to  pick 
various  streamlines  in  the  flow  field  for  the  purpose  of  determining  the 
effect  of  rate  chemistry  upon  electron  density.  Prescribed  pressure 
distributions  along  each  streamline  were  given  as  inputs  to  this  one 
dimensional  calculation,  the  insensitivity  of  the  rate  processes  to 
pressure  distribution  having  been  established  in  Reference  3.  The  gas 
was  assumed  to  consist  of  a  basic  mass  fraction  composition  of  0,  234 
oxygen,  0.  766  nitrogon  and  to  dissociate  into  a  mixture  involving  seven 
chemical  species,  namelyj  N,  N a,  O,  Oa,  NO,  NO+  ,  o“,  which  partake 
in  the  following  reactions: 

Na  +  X  «  2N+  X 
Oa  +  X  *  20+  X 
N|  +  Oa«*  2NO 
NO+X  *  N  +  O  +  X 

NO  +  X  2  NG+  +  e"  +  X 
N  +  0  **  NO+  +  e" 

Na  +  O  *  NO  +  N 
NO  +  O  ~  Oa  4  N 
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Each  individual  component  behaves  as  a  thermally  perfect  gas,  individually 
in  thermal  equilibrium,  although  it  may  partake  in  dissociative  reactions. 
The  details  of  this  analysis  may  be  found  in  Reference  4. 

This  non-equilibrium  molecular  dissociation  calculation  for  the  body 
streamline  was  continued  to  approximately  30  nose  radii  downstream  where 
the  pressure  nearly  approached  the  free  stream  value.  The  results  served 
as  input  to  the  far  wake  analysis  which  is  briefly  discussed  below. 

The  far  wake  analysis  was  carried  out  by  utilizing  the  methods 
given  in  Reference  5  which  involves  an  integral  approach  for  the  solution 
of  the  properties  of  a  laminar,  axisymmetric,  hypersonic  wake  with  coupled 
diffusion  and  rate  chemistry.  This  integral  technique  involves  the  solution 
of  the  applicable  conservation  equations  (mass,  momentum,  energy,  and 
species)  on  the  average  over  a  normal  cross-section,  and  exactly  on  the 
axis. 

A  discussion  of  the  results  of  the  preceding  analysos  may  be  found 
in  the  next  section. 

II.  DISCUSSION  OF  RESULTS 

The  body  streamline  and  five  external  streamlines  wore  chosen  to 
illustrate  the  effects  of  rate  chemistry.  The  relative  location  and  non- 
dimensional  distance  along  each  streamline  may  be  found  in  Figure  2. 
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Table  I  below  lists  the  shock  inclination  and  entropy  for  each  streamline. 

The  geometry  is  explained  in  Sketch  1. 

T^ble  I 

Streamline  Properties 

Streamline  «  rads  Entropy,  S/R 

Upstream  of  Downstream 

Wake  Shock  of  Wake  Shock 


Body 

0i  =  .93046 

47.  490 

47. 856 

A 

1.  1838 

,46.  652 

47. 043 

B 

1.  0691 

45.  000 

45.  364 

C 

1.  0056 

44. 348 

44.  636 

D 

.81319 

41.  793 

41.  796 

E 

. 51777 

37.  690 

m  m  m 

Sketch  1 


-5- 


Figure  3  shows  the  electron  density  distribution  along  the  pellet  body 
streamline,  where  the  initial  specie  mass  fractions  were  taken  at 
their  equilibrium  dissociation  values.  No  comparison  was  made 
with  a  flow  model  which  assumes  that  the  initial  species  concentration 
is  equal  to  the  free  stream  composition  since  previous  studies  have 
shown  the  flow  along  the  body  streamline  to  be  in  local  equilibrium  up 
to  the  vicinity  of  the  sonic  line.  The  isolated  points  indicate  the  local 
equilibrium  values,  1.  e. ,  the  electron  density  which  would  exist  if  the 
ionization  were  in  equilibrium  with  the  local  frozen  concentrations  of 
atomic  oxygen  and  nitrogen.  Note  that  for  small  values  of  S,  the 
electron  density  is  close  to  the  local  equilibrium  value,  but  the  ex¬ 
pansion  corner  causes  considerable  deviation  which  persists  along  the 
remainder  of  the  body  streamline.  The  second  discontinuity  in  electron 
density  occurs  at  the  base  of  the  pellet  where  the  flow  expands  to  the  base 
pressure.  The  third  discontinuity  represents  the  pressure  jump  which 
occurs  at  the  beginning  of  the  wake  shock. 

Figures  4,  5,  6,  7  and  8,  give  the  electron  density  distribution 
resulting  from  the  rate  chemlsty  calculations  along  external  streamlines 
in  th«  flow  field.  Initial  points  for  the  non-equilibrium  dissociation  calcu¬ 
lation  and  the  relative  position  of  each  streamline  may  be  found  in  Figure  2. 
The  dashed  curve  represents  the  electron  density  distribution  resulting 
from  the  initial  species  concentration  taken  at  its  equilibrium  dissociation 
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value.  The  solid  curve  shows  the  distribution  resulting  from  the  initial 
species  concentration  taken  equal  to  the  free  stream  value  (  Qq  =  .  234, 

=  ,7  66)  .  Note  that  for  streamlines  A  and  B  (Figures  4  and  5)  the 
electron  density  distribution  given  by  the  initial  frozen  free  stream 
composition  assumption  rapidly  approaches  local  equilibrium,  and  in 
fact  for  Streaniline  A  overshoots,  and  then  closely  follows  the  distri¬ 
bution  given  by  the  equilibrium  dissociation  assumption,  although 
neither  distribution  can  maintain  the  local  equilibrium  value  for  non- 
dimensional  distances  greater  than  S*»  0,2. 

Quite  a  different  situation  exists  for  streamlines  D  and  E, 
(Figures  7  and  8)  .  Here  the  distribution  corresponding  to  initial  free 
stream  composition  does  not  approach  local  equilibrium  and  is  several 
orders  of  magnitude  apart  from  the  curve  giving  the  electron  density 
resulting  from  the  equilibrium  dissociation  assumption. 

The  results  of  the  nonequilibrium  calculation  for  streamline  C 
represents  a  borderline  case.  For  small  S  the  electron  density 
distribution  does  not  attain  local  equilibrium  although  it  closoly  approaches 
it  and  beyond  S»  0.4  both  distributions  do  not  differ  by  more  than,  a  factor 


of  3. 
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In  summary  then,  the  results  of  the  rate  chemistry  calculations 
point  to  the  existence  of  two  regions  in  the  flow  field  even  for  the  small 
body  geometry  considered  here.  The  first  region  is  characterized  by 
relaxation  distances  which  are  small  when  compared  to  the  nose  radius 
and  includes  the  streamlines  crossing  a  strong  bow  shock.  The  flow 
in  the  second  region  is  characterized  by  large  streamwise  pressure 
variations  over  the  relaxation  distance.  The  flow  entrained  within  the 
first  region  approaches  local  equilibrium  and  atomic  recombination 
freezes  while  the  recombination  of  ions  and  electrons  persists  due  to 
its  faster  rate.  Within  the  second  region  the  flow  never  attains  local 
equilibrium)  the  peak  values  of  dissociation  and  ionization  being  realized 
before  the  onset  of  freeze  due  to  the  interaction  between  the  driving  con¬ 
ditions  for  dissociation  and  the  damping  effect  of  the  rapid  streamwise 
pressure/temperature  decay. 

Figures  9,  10,  and  11  give  the'electron  density  profiles  in  the 
shock  layer  at  the  1st  expansion  corner,  midway  back  on  the  pellet 
( S«»  0.35118)  and  at  the  pellet  base,  respectively.  These  figures  illu¬ 
strate  the  normal  variations  in  electron  density  which  exist  duo  to  rate 
chemistry  considerations. 

Figures  12,  13,  14,  and  15  show  the  axial  distribution  of  electrons, 
temperature,  density  and  velocity  in  the  far  wake  respectively.  The  solid 
curve  represents  the  results  of  the  present  detailed  analysis;  the  dashed 
curve  taken  from  Reference  6  and  included  for  comparison,  shows  the  results 
of  the  approximate  analysis.  While  differences  in  temperature,  density  and 
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velocity  exist  due  to  pressure  differences,  there  is  remarkably  good 
agreement  (Figure  12)  in  the  electron  density  distribution  which  cor¬ 
roborates  the  approximate  analysis  scheme. 

Results  of  the  detailed  rate  chemistry  analysis  of  the  flow 
along  streamlines  in  the  shock  layer  and  near  wake  flows  indicates 
that  due  to  the  small  size  of  the  pellet  the  flow  is  essentially  frozen 
down  to  the  point  where  the  far  wake  begins.  Therefore  since  the  wake 
imput  conditions  are  determined  by  the  body  otreamlino  rate  chemistry 
one  should  expect  only  minor  differences  betwoen  the  initial  conditions 
for  the  far  wake  calculations.  In  the  far  wake,  as  discuoacd  in  Itef.  6, 
due  to  the  small  dimensions  involved,  diffusion  effects  overwhelm  the 
effects  of  rate  chemistry  to  such  an  extent  that  the  peaks  in  electron 
density  predicted  for  a  much  larger  vehicle  and  for  which  rewultu  are 
presented  in  Riferonce  7  are  not  observed. 
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III.  SUMMARY 

The  detailed  supersonic  flow  field  about  an  MIT  firing  range 
pellet  was  computed  by  the  method  of  characteristics.  The  a  priori 
assumption  of  frozen  air  composition  was  made  in  light  of  the  results 
given  in  Reference  3.  This  characteristics  solution  yielded  pressure 
distributions  which  were  used  as  input  to  a  non-equilibrium  molecular 
dissociation  analysis  along  the  body  streamline  and  several  other 
external  streamlines  in  the  flow  field. 

The  effect  of  various  flow  model  assumptions  on  electron 
density  along  each  streamline  is  presented  and  electron  density  profiles 
acToss  the  shock  layer  are  given  for  three  axial  stations. 

The  far  wake  analysis  of  Reference  5  was  used  to  determine 
the  axial  electron  density  distribution  a  few  hundred  nose  radii  downstream 
of  the  pellet,  and  a  comparison  of  the  results  thus  obtained  was  made  with 
those  given  by  the  approximate  calculation  of  Reference  6.  The  agreement 
between  the  results  of  the  actual  and  approximate  calculations  is  excellent. 
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Fig.  L.  Aeroballistic  Pellet  Geometry 
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